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INTRODUCTION 


Natural  electromagnetic  (EM)  radiation  produced  from  surrounding  environments 
and  radio  frequencies  emitted  from  man-made  avionics  sources  interrupt  the  command  and 
communication  signals  necessary  for  desired  combatant  devices.  Knowledge  due  to 
meticulous,  persistent  research  has  produced  state  of  the  art  construction  of  aircraft 
component  systems  with  advanced  composite  materials.  The  achievement  of  structural  and 
economical  goals,  such  as  weight  and  cost,  and  electromagnetic  compatibihty  are  sought 
through  such  theoretical  basis  as  cascading  matrices  and  magiietoelectric  coupling.  Hence, 
a  theoretical  investigation  of  advanced  composite  behavior  in  response  to  various  radio 
frequencies  will  certainly  lead  to  superior  performance. 

The  various  types  of  composite  materials  under  investigation  are:  1)  graphite/epoxy 
composites;  2)  artificially  manufactured  chiral  composites;  3)  quasioptic  composites;  4) 
omega  particulate  composites;  and  5)  dielectric  composites.  Graphite/epoxy  composites 
and  chiral  composites  will  dominate  the  focus  of  this  review  due  to  the  limited  open 
literature  discussion  of  the  other  three  materials  classes. 

In  order  to  accurately  approximate  the  controlling  properties  that  effect  the 
electromagnetic  response  of  composite  materials,  tests  must  be  performed  to  determine  the 
values  of  the  electromagnetic  parameters.  Although  the  various  tests  are  not  discussed,  the 
controlling  electromagnetic  parameters  of  interest  for  graphite/epoxy  composites  are  the 
dielectric  permittivity,  magnetic  permeability  and  resistance  (or  conductivity).  Once  the 
essential  electromagnetic  parameters  values  are  known,  effective  computational  analysis  is 
possible. 

The  controlling  electromagnetic  parameters  for  chiral  composites  are  the  effective 
permittivity,  effective  permeability  and  the  effective  chirality  parameter.  The  study  of  chiral 
composites  in  relation  to  electromagnetic  waves  has  been  relatively  new  due  to  naturally 
occurring  chirality  observed  during  optical  activity.  Chiral,  “left”  or  “right”  handedness, 
refers  to  the  natural  response  of  a  microscopic  particle  to  rotate  polarized  light.  Since  light 
waves  are  electromagnetic  waves,  precise  structural  and  electromagnetic  analysis  can  be 
used  to  predict  of  the  phase  change  induced  by  chiral  composites. 

This  paper  contains  four  major  sections.  Section  I  contains  the  “Foundational 
Theory”  which  structures  theoretical  arguments  such  as  cascading  matrices  via  Maxwell’s 
equations.  Section  II  describes  two  Aeoretical  models  used  to  analyze  graphite/epoxy 
composites:  1)  the  Equivalent-Transmission-Circuit-Line  Model  and  2)  the  Filament- 
Current  Phase-Correction  Model.  Section  HI  provides  a  brief  overview  of  theories  that 
define  the  research  of  helical  and  omega  shaped  particles  embedded  within  a  dielectric 
matrix. 
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1.  Foundational  Theory:  The  Reflection  and  Transmission  of  EM  Waves 

Throughout  the  investigation  of  improving  advanced  composite  materials,  one 
underlining  theory  remains  the  same  -  the  multiplication  of  cascading  matrices.  Matrices 
are  the  three-dimensional,  mathematical  representation  of  the  composite  structure.  A 
“transmission”  matrix  multiplied  by  a  vector  matrix,  composed  of  incident  and  reflection 
coefficients,  yield  a  third  matrix  that  produce  an  increase  or  decrease  in  reflection  and 
incident  vector  coefficients  on  a  fourth  or  following  matrix.  Thus,  composites  are  best 
modeled  as  the  multiplication  of  infinite  plane  matrices  in  the  x  and  y  direction  with 
specified  boundary  conditions  in  the  z  direction.  In  the  following  examples,  the 
mathematical  processes  are  outlined  for  a  single-layer,  infinite  composite  slab  illuminated 
by  a  sinusoid^  electromagnetic  wave  for  normal  and  oblique  incidence. 


Special  Case:  Normal  Incidence/Homogeneous  Isotropic  [12],  [26] 

Region  a  T  x  Region  1  Region  b 


z  =  0  z  =  d 

Figure  1-1 

Plane  wave  equations  representing  the  E-fields: 

Region  a:  E/  =  A/  e  x  and  E/  =  A/  e  x 

Region  b:  Eb'  =  e  x  and  Ej  =  Ab*  x 

Composite:  E/  =  A, '  e  x  and  E/  =  A/  x 

(Region  1) 

Plane  wave  equations  representing  the  H-fields: 

Region  a:  =  1  A^'  e  y  and  H/  =  -1  AJ  e  y 

fio  Tlo 

Region  b:  Hb'=  1  Ab'  e  y  and  Hb'  =  -i  Ab'  e  y 

,  'Ho  .  .  fio 

Composite:  H/  =  1  Aj  *  e  y  and  H/  =  -1  A/  y 
(Region  1)  r[^  Tij 
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An  individual  homogeneous  conducting  layer  can  be  characterized  in  terms  of  the  tangential 
electric  and  magnetic  fields  in  Region  a  and  Region  b.  These  conditions  yield  the  matrix 
equation: 


Ta,'! 

[Ti  “  Ti  '"1 

rvi 

J 

Lt^21  ^^22J 

La;J 

=  g  +  [e  - 

(1 

411, 

(1  +Tli)" 

=  U-Jlif  (e  - 

eik.d) 

4Tli 

'X'  2 1  __  'y  1 2 

Ij  —  -  ll 

Ti ""  =  (I  +  r\i?  [e  -  (1 

4T1i  (1  +  Til)" 


with 

Thus,  the  cascading  of  n  isotropic  composite  layers  between  the  z  =  0  and  z  =  d  interfaces 
are  given  by: 


r^'i 

n  [Tj" 

Ti  1  [a;! 

=  (n 

T,  “  J  La;J 

La,‘ J 

i=lLTi'' 

and  the  ith  layer  transmission  parameters  T^  are  given  by 


T  "  =  a  +Tii)^  [e  - 

4Tli 

(1  +  Tli)" 

T  =  fl  -Tii')^  fe  - 

e^"‘) 

4Tli 

ry  21  _  ly  1 2 

T  =  n  +Ii\.f  [e  - 

(1  -Tl,)"e 

4Tli 

(1  + 

with  Tii=rii/ri„, 
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Special  Case: 


Oblique  Incidence/Homogeneous  Isotropic  [12],  [26] 


Region  a 


Region  1  Region  b 


For  the  oblique  incidence  case,  two  separate  subcases  are  examined,  the  transverse  electric  (TE) 
case  and  the  transverse  magnetic  (TM)  case.  The  equations  that  represent  the  electric  and  magnetic 
field  for  each  case  are  shown  below,  following  are  the  respective  transmission  parameters. 

Transverse  Electric  (TE)  Case 

The  respective  electric  and  magnetic  fields  for  Region  a: 

e;  =  a;  e  -j‘'.(-“9o+xsin0„)  ^ 

E/  =  -a;  e  [X  cosGo  +  z  sinGo] 

Q  -jk„6sm0„+  zcos0„)  ^ 

f|o 

|.]^r  1  A  ’’  e  cos0„  +  X  sm8„)  ^ 

flo  " 

The  respective  electric  and  magnetic  fields  for  Region  b: 

Eb'  =  -A,'  e  [X  00800+  z  sinGo] 

=  a;  e  [X  cosGo  -  z  sinGo] 

_  2  g  -jls.(xsin8„-  (z-d)cos8„)  ^ 

Tic 

J  g  -jko(xsm0„+  (z-d)  cos0„)  ^ 

Tlo 
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The  respective  equations  that  relate  Region  a  with  Region  1  at  the  z  =  0  interface  are: 


B,  =  1/2  (%  +  1/v,)  V  +  1/2  (Ti,  -  1/Vi) 

Q  =  1/2  (Til  -  i/vj  a;  +  1/2  (Ti,  +  1/Vi)  a; 

The  respective  equations  that  relate  Region  1  with  Region  b  at  the  z  =  d  interface  are: 

Bi  =  1/2  (T|i  -  1/Vi)  V  e  +  1/2  (Til  +  1/Vi)  K  e 

Cl  =  1/2  (rii  +  lA^i)  A,'  e  +  1/2  (t1i  -  lA^i)  Ai,^  e 


Therefore,  transmission  parameters,  Tj ,  for  the  transverse  electric  (TE)  case  of  a  composite 
slab  of  ith  layers  are  given  by: 

T  “  =  (I  +  y].y.f  (e  -  a_LUsYi)"  e 

4TliVi  (1  +  TliVi)" 

rj.  12  _  1  _  (e  _  g  jkidiCosft,^ 

4riiVi 

•y  21  _  'T>  12 

1  i 

Ti"'  =  (I  +r|,v,l^  (e  -  (I  -  e 

4TliVi  (1+riiVi)" 


where  0i  =  sin  ’  (ko/ki  sin  0o) 

Tli  =  Tli/T|„ 

Vj  =  COS0i  /  COS0O 
Transverse  Magnetic  (TM)  Case 
The  respective  electric  and  magnetic  fields  for  Region  a: 


i  e>.  +xsin0„) 


g  r  _  g  -jk„(-z  COS0„  + X  sin9J  , 


H,'  =  1  A,'  e  [-X  COS0O  +  z  sin0o] 

=  1  a;  e  cos0„  +  xsin0„))^^  ^  2 

Tlo 
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The  respective  electric  and  magnetic  fields  for  Region  b: 

gj  _  g  -jk„(-(z-d)  cose„  +xsinej  y 

Q  -jk„((z-d)  sine„+xsinej  y 

"Ho 

=  1  e  cos0„.xsmej|-_^  ^  ^ 

Tie 

The  respective  equations  that  relate  Region  a  with  Region  1  at  the  z  =  0  interface: 


B,  =  1/2  (1  + 11,  /V,)  a;  +  1/2  (1  - 11,  /V,)  a; 

C,  =  1/2  (1  - 11, /V,)  A.‘  +  1/2  (1 +11, /v,)  A.' 

The  respective  equations  that  relate  Region  1  with  Region  b  at  the  z  =  d  interface: 


Bi  =  1/2  (1  -  Til  /Vi)  A,'  e  +  1/2  (1  +  Til  A'l)  e 

Cl  =  1/2  (Til  /Vi  +1)  A  e  +  1/2  (-  Til  /Vi  +1)  A'  e 

Therefore,  transmission  parameters,  T^ ,  for  the  transverse  magnetic  (TM)  case  of  a 
composite  slab  of  ith  layers  are  given  by 

T  “  =  fri,  +  (e  -  fri,  -  e 

4TljVi  (Tli  +  V.f 

=  (tIi  -  y.l^  (e  _  e  -jkid,cosft,^ 


4TliVi 


12 


fri,  +  Vjl"  -  rn,  - 

4TliVi  (Tli  +  Vi)" 


where  0i  =  sin ''  (koA  sin  Go) 

n,  =  r[J^o 

l/Vj  =  cosGo  /  cosGi. 

Conesponding  formulas  for  the  anisotropic/  normal  incidence  and  anisotropic  oblique 
incidence  are  given  in  reference  [26]. 
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11.  Graphite/Epoxy  Composites 

The  first  composite  class  discussed  in  this  paper  are  anisotropic  carbon  fiber 
composites  (CFC).  T^e  carbon  fiber  composites  (Figure  1-2)  normally  utilized  in  most 
aerospace  applications  are  the  graphite/  epoxy  (G/E)  composites.  The  special  properties  in 
specific  directions  are  designed  by  choosing  the  required  fiber  and  matrix.  However,  the 
most  important  mechanical  properties  are  the  specific  tensile  strength  or  specific  modulus. 
These  properties  are  the  central  to  design  decisions  for  aircraft  applications  because  they 
present  the  ratio  of  strength  (or  stiffness)  to  weight  [12]. 


Figure  1-2  Schematic  Diagram  of  CFC  [12] 


In  addition  to  the  composite  fiber,  the  composite  matrix  properties,  in  the  form  of 
resin,  performs  two  important  tasks:  (1)  separation  of  fibers  to  maximize  their 

contributions  to  strength  and  minimize  the  possibility  of  cracks  and  buckling  caused  from 
fibers  in  contact  and  (2)  enabling  the  transfer  of  stresses  between  fibers  by  sufficiently 
binding  the  fibers  throughout  the  material. 


The  following  two  sections  outline  the  basic  formulation  of  the  Equivalent- 
Transmission-Circuit-Line  (ETCL)  model  and  the  Filament-Current  Phase  Correction 
(FCPC)  model  for  graphite/epoxy  composites. 


II.  A.  The  Equivalent-Transmission-Circuit-Line  Model 

An  equivalent-transmission-circuit-line  model  represents  the  electric-  and  magnetic- 
field  responses  of  graphite/  epoxy  composites  to  radio  fequencies  as  a  complex  yet  simple 
circuit  of  alternating  current.  In  addition  to  the  electromagnetic  properties  of  permittivity 

(£  ),  permeability  (/i  )  and  conductivity  (cr  );  inductance  (L  ),  capacitance  (C  )  and 
conductance  (G  )  properties  are  utilized  to  describe  transient  propagation  of  plane  waves 
through  composite  structures.  Hence,  graphite/  epoxy  composites  are  characterized  as  a 
combination  of  semiconducting  and/or  insulating  material. 
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The  wave-transmission-matrix  (WTM)  method  was  initially  used  to  construct  the 
solution  of  frequency-domain  Maxwell’s  equations  for  laminated  anisotropic  composites 
illuminated  by  an  obliquely  incident  plane  wave.  By  imposing  the  necessary  boundary 
conditions  on  the  interfaces,  and  finally  cascading  the  wave-transmission  matrix  for  each 
lamina,  a  derivation  was  developed  for  the  relationships  between  the  incident,  reflected, 
and  transmitted  field  [24],  In  the  following  problem  statements,  a  brief  analysis  of 
anisotropic  homogeneous  and  anisotropic  inhomogeneous  laminates  are  given  for 
graphite/epoxy  composites. 

11.  A.  1.  Frequency-Domain  Analysis  of  Graphite/Epoxy  Composites 
11.  A.  1.  a.  Special  Case:  Anisotropic  Homogeneous  Laminates  [24] 
Problem  Statement 

Consider  a  composite  structure  with  M  number  of  layers,  each  described  by  discrete 
electromagnetic  properties  such  as  permeability  ,anisotropic  permittivity  e  ,  and 
conductivity  O (Figure  1 ,  Appendix  A).  Where  the  permeability  is  essentially  uniform, 
while  the  permittivity  and  conductivity  are  expressed  as  tensors  in  the  following  form: 


ol 

LO  0  aj„ 

UC/O  Xj/  -r  C-y  JfcAt  Xj/ 

£yx  =  (^x-  ~  £y)  cos  $sin  0 
siri  <p  +  £y  coi  0 

e  =  £. 

zz  z 

The  primed  parameters  of  permittivity  are  the  composite  principal  axes  and  (])  is  the  angle 
between  the  global  and  the  principal  axis  (Figure  2-2).  The  expressions  for  the 

conductivities  are  similar  to  those  used  for  the  respective  permittivities  (cr  ,  <7  =  a  ,  <7  , 


0  1 

~  ^yx  ^yy  ^ 

Lo  0  £j, 


Here 
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Figure  2-2 

Global  coordinates  (black 
coordinate  axis)  and  principal 
coordinates  (gray  coordinate  axis) 
for  laminated  composites  [24]. 


The  circuit  simulation  program,  PSPICE,  has  been  used  to  perform  the  numerical 
analysis  of  the  ETCL  model.  However,  before  subjecting  the  ETCL  model  to  a  simulation 
program,  the  transmission  lines  must  be  individualized  and  their  values  estimated  by 
cascading  3t-  and  T-circuit  sections  as  either  jt-T-jt  or  T-Jt-T  circuit  models  (Figure  2, 
Appendix  A).  According  to  reference  [24],  the  jt-T-jr  circuit  model  produces  better 
convergence  and  numerical  stability  than  the  T-Jt-T  circuit  model. 

In  the  analysis  of  the  ETCL  model  a  composite  lamina  of  length  d  contains  an  odd 
number,  2m  -f  1,  of  jt-  and  T-  circuit  sections,  where  m  is  the  number  of  laminae  of  length 
d  within  the  single-composite  layer.  Assuming  the  same  method  applied  to  inhomogeneous 
composites,  an  accurate  convergence  of  the  n-T-n  circuit  model  may  be  obtained  by 
calculating  N  number  of  circuit  sections  (where  N  =  2m  +  1 .  The  total  length  of  the  n-T-n 
circuit  section  is  characterized  by  a  length  1  where  the  Jt-section  (1,)  is  half  the  length  of  T- 
section  (Ij)  of  the  circuit  [15].  For  the  ETCL  model,  the  length  of  the  jt-  and  T-sections 
should  less  then  half  the  effective  skin  depth,  where  the  skin  depth  is  a  function  of 
effective  conductivity  and  frequency  [24]. 

=  V  _ 1 _ . 

The  appropriate  coupled  transmission  line  equations  for  the  homogeneous 
anisotropic  ETCL  model  are: 


6  0)  =  -L  •5_t_ 

6z  5t 

5_.-  =-(C  L+G  )  •  0) 

6z  6t 

The  equivalent  inductance,  capacitance,  and  conductance  matrices  in  the  coupled 
transmission  line  equations  for  the  homogeneous  case  are: 

fl  -  sin^9  01  \e^  1 

C  =e^ 

L  0  1  J  l-£^  Syf-  sin^e  J 
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the  necessary  boundary  conditions  for  the  voltage  and  current  are  given  by; 


)  +  77 p  ip  (0,t  )  =  2\)p  ‘(t  )  at  incident  boundary 

where  p  =  1 1  or  ± 

t  )  +  r]  p  ip  (D,  t  )  =  0  at  transmitted  boundary 

Where  p  denotes  the  parallel  component  or  the  perpendicular  E-field  of  the  plane  wave  to 
the  x-z  plane  of  incidence.  The  voltage  and  current  wave  vectors  are  defined  as  in  [24]: 


\)(  z,  O  =  Td  1 11  =  r  5^  1  =  J  r  £,  (  z;  fy  )  do) 

LvJ  1-BX^o  ^  1-E^{2;co)\ 

c{z,t  =  =J  \h^{z;co)  ^d^dco 


LiJ  i-^J 


y-^x=0 


L  -//,  ( z;  ft)  )  J 


Here  the  transient  waveform  as  a  function  of  the  incident  angle  0  for  the  plane  wave  can  be 
expressed  as: 


S(x,  z,  t  )=1  E  (z;o)  )e’‘^ 

1^{x,z,t)=-]  H  (z;  CO)  dco 

where  k„  =  (oV|i„e  „  =  (O  /  c  is  the  free-space  wavenumber  and  c  is  the  velocity  of  light. 

The  combination  of  the  ETQ^  model  and  the  wave-transmission-matrix  (WTM) 
method  can  be  used  to  predict  the  transmittivity  and  reflectivity  terms  which  are  basically 
mathematical  expressions  of  the  transient  plane  wave  behavior  within  a  composite 
structure.  In  the  PSPICE  simulation  program,  the  command  “.AC”,  conducts  the 

frequency-domain  analysis.  The  transmittivity  coefficients  T||  j|,  T||  Tj,  j^,  and  Tj^jj  and 

the  reflectivity  coefficients  R|  |  j  | ,  R|  |  Rx,x’  R^ j  |  in  the  frequency-domain  are  given 

by  the  following  definitions  [24]: 


T 


P.‘! 


An  estimation  of  the  copolarization  term  for  transmittivity  in  the  lower  frequency  range  for 
specific  case  where  each  lamina  has  equal  thickness  is  given  by: 


T 


p.p 


2 


where  77^  ,  the  characteristic  impedance  of  the  incident  wave,  is  given  by  77 1 1  =  77  ^  cos0 

and  77  j^=  77  ^  /  COS0  for  the  E|  |  and  Ej^  wave,  respectively;  and  ti  ^  „  is  the  intrinsic 

impedance  of  the  free  space  [24]. 


a 


eflf 


M 

J-TZ  (a^)J 
M  J 
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(  (ohm/m)  is  the  effective  conductivity  of  the  mth  ply  within  the  composite  structure. 

This  can  be  described  by  the  corresponding  formulas  for  the  E|  |  and  wave  respectively 
as: 


cos^  (])  +  (Oy.)^  sin^  ([) ,  for  E, ,  wave 

(  ^e£F  )m  “ 

(a^.)„  sin^  (j)  +  cos^  ^ ,  for  E^  wave. 

Note:  Although  the  formula  for  the  cross-polarization  terms  are  not  given,  when  a  wave 
performs  multiple  reflections  in  a  region,  the  cross-polarization  terms  will  be  dominant  if 
the  thickness  of  the  region  is  near  a  quarter- wavelength  [24]. 

11.  A.  1.  b.  Special  Case:  Anisotropic  Inhomogeneous  Laminates  [15] 

Problem  Statement 


Consider  a  composite  structure  with  one-layer  of  length,  D.  The  electromagnetic 
properties  of  the  slab,  such  as  permittivity,  permeability  and  conductivity  are  described  as 
tensors  and  expressed  in  the  following  form: 


£  = 


[e  £  0  1 

'  xc  xy  ' 

^yy  ^ 

L  0  0  ej 


I^yy  0 

LO  0  eJ 


o1 

^yx 

0 

Lo 

0 

o-J 

These  tensors  are  assumed  to  be  a  function  of  z,  but  independent  of  frequency.  The 
governing  equations  for  the  transverse  field  components  are  given  by: 


5  t)  =  -L  •  5  c 
5z  5t 

5_c  =-(C  5.+  G  )  ■  -u 
5z  5t 


Where 

ri 

-  sin^G 

o1 

-^,1 

II 

fizz 

C  =£„ 

L 

0 

1  J 

sin^e  J 
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G  = 

L-a^  o^J 

The  incident  plane  is  chosen  to  be  the  x-z  plane  with  an  incident  angle  of  0  and  arbitrary 
transient  waveform.  The  plane  wave  is  divided  into  two  components  which  are  the  parallel 

and  perpendicular  electric-field  in  reference  to  the  plane  of  incidence,  represented  by  E|  |  and 
respectively.  Considering  the  fields  at  the  x=0  plane,  the  voltage  and  current  wave 
vectors  are  defined  as: 


v(z,  t  )  =  Fi)  |)1  =  {  £^~\ 

L'Oj.J 

i{z,  t  )  =  =  \ 

LiJ 

Last  and  most  important,  the  necessary  boundary  conditions  can  be  stated  as  follows 
)  +  ty  p  ip  (0,t  )  =  2X)p  ‘(t  )  at  incident  boundary 

where  p  =  1 1  or  ± 

Vp(D,t  )  +  ri  p  Cp  (D,t  )  =  0  at  transmitted  boundary. 

11.  A.  2.  Time-Domain  Analysis  of  Graphite/  Epoxy  Composites 

Special  Case:  Incident  Transient  Pulse,  EMP 

A  specific  type  of  electromagnetic  radiation  produced  from  surrounding  natural  or 
man-made  environments  is  the  electromagnetic  pulse  (EMP).  For  this  case,  consider  the 
field  B{z,  t  )  at  x=0  plane,  under  an  obliquely  incident,  electromagnetic  pulse  waveform 
S  ^(t )  given  by 


S\t  )  =  E^{e-^-e-P')u{t  ) 

here  u{t  )  =  1  for  r  >0.  Also,  the  values  of  =  l/(e““'°  -  )  and  =  In  (Pla)/(P-a) 

[24].  When  plotting  incident  electromagnetic  pulse  versus  time  or  frequency,  the  values  of 

a  and  P  determine  the  shape  and  frequency  spectrum  of  the  EM  pulse. 
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11.  B.  The  Filament-Current  Phase-Correction  Model  [14] 

The  goal  of  the  filament-current  phase  correction  model  is  to  reduce  the 
unpredictable  scattering  patterns  which  disrupt  communication  and  control  signals  utilized 
by  various  applications.  The  electric-field  radiated  by  a  fiber  embedded  within  the  dielectric 
matrix  and  the  electric-field  resulting  from  the  incident  plane  wave,  constitute  the  total 
electric-field  produced  by  the  composite  structure.  The  electric-field  resulting  from  the 
incident  plane  wave  can  be  expressed  as: 

£  =  exp  (-ikx  sin  6  +  i  ky  cos  0) 

while  the  scattered  electric-field  resulting  from  the  fiber  embedded  within  the  dielectric 
matrix  can  be  expressed  by  the  equation 

=  - krjlj.  D(x,y) 

4 


where  T]  =  '4nJ(£^£^ )  and 


D  (x,  y)  =  Y,  2  i  exp  ( -  a^JyJJ  exp  (  -  ik„x ) 


and 


here 


k^  =  ksin  6  +2nK 
d 


= 

^/k^^  -H,  if  <1^^  . 

In  this  analysis  the  composite  structure  is  described  by  a  periodic  fiber-matrix 
(Figure  3-2).  The  total  electric-field  can  be  substituted  into  an  equation  representative  of 
the  concentrated  electric  current  along  the  center  of  the  conducting  fiber. 


/y  =  II J  dxdy  =  If  oE  dx  dy 

A  A 

where  E  is  the  total  field  and  A  is  the  fiber  area. 
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From  this  approximation  of  the  electric  current,  the  reflected  and  transmitted  electric-field 
can  be  computed. 


oil  exp  (  -  i  k X  sin  6  +  i  ky  cos  6 )  dx  dy 


1  +  akr)  HE 


exp  (  -  i  k^x)  dx  dy 


The  relation  between  the  reflected  and  transmitted  electric-field,  approximated  electric 
current  and  phase-correction  model  follow  in  the  succeeding  discussion. 

The  phase-correction  model  follows  the  logic  of  the  foundational  theory,  the 
multiplication  of  matrices.  In  reference  to  the  equation  for  the  scattered  electric-field,  the 
Fourier  transform,  D(x,  y  ),  contains  an  infinite  number  of  discrete  terms  numbered  as 
individual  Hoquet  modes.  The  computations  for  the  phase-correction  model  are  formulated 
for  a  single  Floquet-layer  within  a  composite  structure  (Figure  4-2),  where  the  Roquet 
layer  contains  2m  +  1  Roquet  modes. 

The  incident  wave  that  illuminates  the  composite  on  the  initial  fi-ee  space/  dielectric 
interface  is  represented  by  the  vector  Aj 

A,  ... 

Where  E,  (I  =  -m,  ...  0,  .  .  .  ,  m)  is  the  amplitude  of  the  /  th  Roquet  mode  and  the 
superscript  r  denotes  the  transpose  [14].  The  same  principle  follows  for  the  vectors  Aj,  . 

.  .  ,  Ag  and  Bj,  .  .  .  ,  within  the  2m  -t- 1  Floquet  modes. 

Thus,  the  reflection  and  transmission  phenomena  at  the  air-matrix  interfaces  may  be 
represented  by 


[BJ^  =  Ta,]  , 

iB,}  Vif,  LaJ 


[b,!  =  Fa,! 

IbJ  It^rJ  LaJ 

and  those  at  the  Floquet  layer  (fiber  grating)  by 


[5,1  =  [R  T^  [a,! 

IbJ  It  r]  LaJ 

The  subscripts  fd  for  the  R  and  T  coefficients  are  the  respective  diagonal  reflection  and 
transmission  matrices  from  free  space  to  the  dielectric  matrix,  and  the  subscripts  df  for  the 
R  and  T  coefficients  are  the  respective  diagonal  reflection  and 
transmission  matrices  from  the  dielectric  matrix  to  free  space. 
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The  relationship  between  the  wave  vectors  A,,  and  5,  (where  i  =  1,  2,  ...  5,  6)  and 
the  phase  shift  may  be  expressed  by  the  multiplication  of  the  following  matrices: 


rp,i  = 

[p  + 

0  IFaJ, 

LaJ 

Lo 

p-JLfiJ 

U,!  = 

Fp- 

0  IFaJ. 

\-bJ 

Lo 

p+JLbJ 

Where 

P  +  =  a  diagonal  matrix  with  dimensions  (2m  +  1)  x  (2m  +1),  whose  iith 
element  is  exp  [t^,^^^y2] 

P  -  =  the  complex  conjugate  of  P  + 

Let  a  combination  of  the  above  matrices  be  represented  by  the  four  full  matrices  M,j, 

Mgj,  andMgg  (analogous  to  the  transmission  coefficients  Tj",  Tj^',  and  T.^^).  The 
combination  of  the  above  matrices  are  represented  by  cascading  the  vectors  [B,  and 
[Aj  expressed  by  the  following  relationship 

fP/l  =  \m,,mJ\a;] 

LpJ  Im,,mJIa,\  [14]. 


Thus,  the  ETCL  model  and  the  FCPC  model  parallel  in  three  different  areas:  1)  Theoretical 
Basis;  2)  Structural  Modeling  and  3)  Fundamental  Compartmentalization.  The  theoretical 
basis  for  both  models  are  the  same,  cascading  matrices.  The  structural  modeling  for  the 
FTCL  model  is  an  AC  circuit,  while  a  periodc  fiber-matrix  describes  the  FCPC  model. 
Finally,  the  fundamental  compartmentalization  for  the  FTCL  model  and  the  FCPC  model  is 
2m  -I-  1  T-  and  k-  circuit  sections  and  the  2m  -i-  1  Floquet  modes,  respectively.  This 
information  has  been  tabulated  in  Table  4. 


Table  4  Summary  of  ETCL  and  FCPC  model 


ETCL  Model 

FCPC  Model 

Theoretical  Basis 

Cascading  Matrices 

Cascading  Matrices 

Structural  Modeling 

AC  Circuit 

Periodic  Fiber-Matrix 

Fundamental 

Compartmentalization 

2m  -1- 1  T-  and  Jt-  circuits 

2m  +  1  Floquet  modes 
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Phase-Correction  Model  for  Graphite/Epoxy  Composite 


interface  1 


interface  3 


Figure  4-2.  Phase-correction  model  for  multi-mode  reflection  and  transmission  due  to  fiber  grating 
within  a  dielectric  matrix. 


III.  Chiral  Composites 
in.  A.  Brief  History 

In  1812,  Biot  observed  the  phenomenon  of  optical  rotation,  or  more  appropriately 
termed,  optical  activity.  Optical  activity  occurred  when  white  linearly  polarized  light 
changes  to  a  rainbow  of  colors  when  shown  through  the  crystal,  quartz  [21],  In  1848, 
Pasteur  accredited  the  rotation  of  plane  polarized  light  and  the  observation  of  the  two  forms 
of  tartaric  acid  salts  to  the  handedness  of  the  material’s  microstructure  [25].  The  term  chiral 

comes  from  the  Greek  word  xpo,  which  means  “hand”  or  a  typical  mirror-asymmetrical 
object  [9].  In  the  early  1900’s,  Lindman  discovered  the  handedness,  chirality,  of 
artifici^y  constructed  materials.  His  measurements  of  plane  polarization  through  a 
collection  of  wire  helices  in  a  lower  frequency  range  have  established  studies  of  observed 
rotation  and  absorption  of  electromagnetic  waves  in  dielectric  material,  since  light  is  an 
electromagnetic  wave  [21],  [25].  In  1956,  Winkler  performed  similar  experiments  with 
200  randomly  oriented  right-handed  copper  hehces  embedded  in  a  dielectric  matrix,  but 
concluded  the  measurements  of  plane  rotation  were  due  to  diffraction  effects  rather  than 
optical  activity  [10],  [21].  He  found  that  over  a  frequency  band  from  1.5  to  2.7  Ghz  the 
electric  field  of  a  linearly  polarized  wave  was  rotated  by  angles  varying  from  0  to  12 
degrees. 


I^ter  in  1957  a^d  1960,  Tinoco  and  Freeman  conducted  numerous  measurements 
of  effectively  anisotropic  chird  material.  They  discovered  that  illumination  of  a  composite, 
consisting  of  right  handed  helices  oriented  parallel  to  the  direction  of  propagation,  yields  an 
electric  field  with  an  angle  variation  from  0  to  30  degrees  in  the  frequency  range  of  8.3  to 

I I. 5  Ghz.  The  resonance  band  frequency  of  12.5  Ghz  yields  rotated  angles  varying  from 
30  to  180  degrees  [10].  Winkler,  Tinoco,  and  Freeman  aU  assumed  the  light  remained 
linearly  polarized  after  leaving  springs,  since  they  lacked  the  equipment  used  to  measure 
copolarized  and  cross-polarized  components  of  the  transmitted  elliptically  polarized  waves 

[10] .  In  1974,  Bohren  introduced  the  decomposition  of  electromagnetic  waves  into  left  and 
right  circularly  polarized  fields  [21].  Therefore,  the  wave  leaving  a  medium  of  springs,  of 
the  same  handedness,  was  usually  elliptically  polarized  [10]. 

Since  then,  the  discovery  of  chiral  helices  has  expanded  to  inclusion  shapes  such  as 
omega  particles  [17].  Once  a  solid  theoretical  and  economical  basis  has  been  formed  chiral 
comp>osites  might  be  utilized  as  a  advance  in  the  control  of  electromagnetic  interference. 
The  following  section  is  a  brief  synopsis  of  the  various  theories  related  to  chiral 
composites. 

III.  B.  1.  Theoretical  Approaches 

Numerous  approaches  have  been  developed  to  account  for  the  behavior  of  chiral 
composites.  The  first  approach  discussed  in  this  section,  the  Waterman-Truell  approach 

[11]  analyzes  the  dielectric-dielectric  composite  in  which  chiral  particles  are  made  of  chiral 
material  but  do  not  necessarily  possess  chiral  shape.  This  approach  assumes  a 
homogenized  medium  with  no  interaction  between  the  microscopic  particles.  The  argument 
begins  by  stating  equations  for  point  dipoles  at  the  origin  and  then  relating  the  scattered 
electric  and  magnetic  fields  to  the  average  dipole  moments. 

From  previous  studies,  medium  characterized  by  the  above  relations  is  isotropic 
chiral.  As  mentioned  before,  propagation  in  chiral  material  produces  elliptically  left  (+)  and 
right  (-)  circularly  polarized  plane  waves.  The  circularly  polarized  plane  waves  for  the 
respective  incident  electric  and  magnetic  fields  are: 
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=  (1/V2)(u^  ±  iUy)  exp  (ik^z) 

^incW  =  ±(l/iT|V2)(u^  ±  iUy)  exp  (ik^^z) 

therefore  the  average  dipole  moments  for  the  circularly  polarized  plane  waves  are  given  by: 

p  =(l/V2)(a^±a,,AlJ(u,±iUy) 

m  =  ±(l/iTi„V2)(  ±  (u,  ±  iUy). 

By  examining  the  electric  field  in  the  far  zone,  the  following  limit  can  be  stated  as: 

E(r)  =  03{cDpJp  -  uXp  •  u^)]  -  k„u,  x  m}{exp  (ik„r)/47ir} 

The  forward  and  backward  plane-wave  scattering  amphtude  can  then  be  described  bv  the 

ennafinns-  J  ^ 


Fforw  -  (ca'/47u)(l/V2)(u,^  ±  iu,)  {p„  (a^  ±  a,,/nj  8  ± 

Eback  =  («'/47c)(1/V2)(u,  ±  iUy)  {p„  (a^  ±  -  8  .(a,,  ± 

^d  the  wavenumbers  for  the  left  (+)  and  right  (-)  polarized  wave  according  to  the 
Waterman-Truell  approach  are:  * 

Ywt  ±  =  K^{  1  +  (coN/k/)[a)(  p^a^  +  8  ±  2  k^a  J 

+  (CON/  kf[a^  oc,,  -h  a,,  (a^  p,  +  oc,,  /  p J] ) . 

The  Maxwell  Garnett  approach  [11]  derived  constitutive  relations  for  the 
homogemzed  chiral  medium  in  the  form  of  bulk  parameters.  The  effects  of  electromagnetic 
radio  waves  was  examined  in  the  low  frequency  range.  Each  chiral  particle  has  been 
encapsulated  m  a  bianisotropic  sphere  which  is  described  by  the  Tellegen  chirality 
pyameter.  Thus,  isotropic  chiral  medium  adhere  to  the  following  Tellegen  constitutive 
relations  given  by  the  respective  electric  and  magnetic  flux  densities 

^  “  E^g  E  +  H 

E  ~  Mmg  H  E 

where 

Emg  =  E  o  J=^2N!(a^^ 

[9p„  8  „  -  3N(8  -H  2p„a  J  -t-  2N"(a^a,,-  a,,")] 

IS  the  permittivity, 


Pmg  -  Po  XS|io_E-o_=i  3N(p^a^^r_28.^ai^fa)jr_ 


[9Po  Eo  -  3N(8„a,,  +  2p„aJ+  2N^(a^a,,-  a  ,2)] 
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is  the  permeability, 
Cmg  “  ^ 


[9Po  ^  o  “  3N(8  +  2p„a  J  +  2N^(a^ai,h 

k  the  Telleeen  chirality  parameter.  The  wavenumbers  for  the  left  {+)  and  right  (-)  circularly 
p£  wavffi^  medium  accorfing  .0  d,e  MaxweU  Game,,  approach  are. 

YmG  ±  =  Ci[^(£MC^MG)  -  ^MG  !  ^  1 

Therefore,  the  intrinsic  impedance 

^MG  “  ^(PmG^Mg) 

impedance  of  the  surrounding  atmosphere. 

III.  B.  2.  Theoretical  Models,  Formulas  and  Equations 

The  generalized  Maxwell  Garnett  model  [2]  for  isotropic  chiral  composites  yields 
the  foUowinf  for  the  effective  electromagnetic  parameters  of  permittivity,  permeability  and 

chirality. 

^  0+  -3-  ^ee“  —A-J 


D 


SPPo 


Peff  =  M}ic+  _i_  (N  — A_) 


D 


See, 


Kjff  —  - j -  N  _^an- 

D  Ve„ 


Po 


where  D  =  1— N  a^—  N  _a^  + 

3ee„  3pp„  9pp„E£. 


A  ~  ^ee  ^imi 


where  e  and  p  are  the  relative  parameters  of  the  matrix  (e^  -  e/e„  p,.  p/Po  [I])- 


The  constitutive  relations  for  the  bi-anisotropic  chiral  media  in  terms^of  el^tric^md 


magnetic 

equations 


ine  consiliuuvc  iciau^uo  xwu  - - - — r -  tV.f»  fr>11r>winCT 

ic  flux  densities  and  corresponding  dipole  moments  are  expressed  in  th  g 
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D  =£.E-iV(£„jiJ  K^.  H 


B  =p.  H+iV(£„pJK.E 
and 

p  =  •  E  +  •  H 

m=a^«H-a^*E 

In  addition  to  the  above  approach,  another  formula  was  developed  for  the  design 
chiral  composites.  The  modified  Maxwell  Garnett  formula  for  chiral  composites,  the  Chiral 
Maxwell  Garnett  Mixing  formula  [6],  describes  the  macroscopic  electromagnetic 
parameters  of  permittivity,  permeability  and  chirality  with  the  necessaiy  electric  and 
magnetic  polarization  terms.  The  argument  begins  by  stating  the  material  constitutive 
relations  which  represent  magnetoelectric  coupling  for  the  bi-isotropic  media 

D  =££  -  iKV(p„£^/f 

B  =pf/-t-iKV(p„£^£ 


Where 

D  =  electric  flux 

B  =  magnetic  flux 

£  =  permittivity  of  chiral  material 

£  „  =  permittivity  of  the  vacuum  (or  achiral  background) 

p  =  permeability  of  chiral  material 

p^  =  permeability  of  the  vacuum  (or  achiral  background) 

K  =  dimensionless  parameter  of  chirality 

Thus,  the  homogeneity  of  the  composite  material  can  be  defined  by  the  effective  parameters 
of  the  mixture.  The  constimtive  equations  described  by  the  average  flux  densities  and  the 
associated  electric  and  magnetic  fields  are: 

0)  =£,jj£-  iK^V(  p„£o)H 
(B)  =  Peff^  +iKeffV(  p„  E^E 

Including  the  electric  and  magnetic  polarizabilty  of  the  constitutive  relations  yields 
(D)  =£„E  +P,  =  e.E  +P„  +  P^ 

(B)  =  yi^H  +  P„=]i„H  +  P^  +  P^ 


The  definition  of  the  polarizability  can  be  stated  as  follows: 
Pu  =  ^PiJ 
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where  i  and;  are  any  combination  of  e  andm  and  let  n  be  the  number  of  spherical  chiral 
scatters  per  unit  volume,  whose  parameters  are  e,  p,  and  k.  Therefore,  the  dipole 

moments,/?^.  ,are  equal  to  or  directly  proportional  to  the  polarizability  factor,  ,  and  the 
Lorentzian  fields,  and  ,that  encompass  the  surrounding  polarization  of  the  field 
created  by  the  shape  of  the  scatter. 

Pee  =  «„ 

Pen,  = 

Pnnn  = 

Pnu  = 

The  exciting  field  for  spherical  inclusions  are  defined  as: 

E,=E  +P^ 

3eo 

3lio 


Polarizability  factors  for  single  chiral  scatters  defined  in  compliance  with  quasistatic 
analysis  yield 

=  47t  a '  e  „  (ii+  2p„J(£.;^£  p,_£  „ 

(|i  +  2p„)(e  +  2e  J-K^  p„  e„ 

47taN(  p,  ej _ ^j3Kp„_£„ _ 

(li+2pj(e  +  2e  J-K^  p„  e„ 


a_  =  4tc  a  J  ^_^_e  „ 

(p+2pj(e  +  2e  J-ie 

=  V(p„  e  J _ j3Kp„_e„ _ 

(P+2p<,  )(e  +  2ej-ie  p„e„ 


where  a  is  the  radius  of  the  sphere  with  em  parameters,  permittivity,  permeability  and 
clurahty.  Thus,  the  effective  electromagnetic  parameters  of  the  homogeneous  chiral  slab 


22 


eeff=eo  +  3/eoX  -  /  (|i  -..tioILziKi4i.-£  oOiii) 

[  (^i  +  2p„)  -  /  (^i  -  p  J] [  (e  +  2e  J  -  /(e  -  e  „)] 

-k"  p„e„(l-/)^ 


Peff  =  Po  +  3/Po  X 


_ (e  + 

[  (p+  2pJ-/(p-pJ][(e  +  2e  J-/(e-e„)] 

-k"  p„e,(l-/f 


Keff=  - 9iK^^ _ 

[(\i+  2pJ-/(p-p„)][(e  +  2eJ -/(£-£„)] 

p„eo(l-/)^ 


where  f  =  n  4tui  ^/3. 

According  to  reference  [5],  there  are  four  additional  constitutive  equations  which 
describe  the  effective  properties  of  isotropic  chiral  composites:  1)  the  Tellegen  equations 
(also  used  to  support  the  previously  cited  Maxwell  Garnett  Approach);  2)  the  Drude-Bom- 
Fedov  equations;  3)  the  Condon’s  Time-Harmonic  equations;  and  4)  the  Post  equations. 
Assuming  the  volumetric  inclusion  portion  of  the  chiral  material  to  be  represented  by  small, 
mutually  noninteracting  spheres.  LetiV  be  the  number  of  spheres  per  unit  volume.  Also 
assume  the  composite  medium  to  be  a  homogeneous,  nonvacuous  region  with  polarization 
and  magnetization  for  the  respective  polarization  and  magnetization  field  approximated  as 

P=Np„,  M=Am„. 

Therefore  the  D  and  B  fields  for  the  effective  medium  can  be  written  as 


D  =  (£„  +  A  aJE  +  i  H , 


These  constitutive  equations  are  applicable  to  isotropic  chiral  composites  [5]. 


The  Tellegen  Constitutive  Equations 

The  Tellegen  constitutive  relations  describe  the  discrete  random  composite  as  an  effectively 
isotropic,  reciprocal  chiral  medium  [5]. 

D  =  £  '^^E+  H 

B  H-  ^ 

Where 

=£„+  N  a,,  =f^o  +  ^  ^  =  iN  ccem 

The  Drude-Born-Fedorov  Constitutive  Equations 

The  Drude-Born-Fedorov  constitutive  equations  are  yet  another  representation  of  the 
effective  parameters  of  isotropic  chiral  medium. 

B  =  pi  xH] 

where 

=  (£,  +  W  a„)  -  (N  a,S I  +  Na^), 

M  +  N  a^)-(.N  ay/ie,+  N  o„), 

=  (N  aJo»/  l(e,+  N  aj(u.  +  N  a^)-  {N  aj=] 

If  the  small  spheres  are  not  chiral,  A  P  are  all  equal  to  zero,  then 

^(DBP)  =  =  e  j  1  +  4jt  Na  -  1)/  (e^  +  2)] 

pi  =pi,{\+AnNa  ^(pi,  -  1)/  (m,+  2)] 

and  the  effective  medium  is  also  achiral  [5], 

The  Condon’s  Time-Harmonic  Constitutive  Equations  [5] 

D  E  +  H 

B  =  pi  E 

where 

e<^^=E^+N  =  N  ajco 


24 


The  Post’s  Constitutive  Equations  [5] 

D  +  B 


H  =  (ll^  B  +  E 

where 


=  ie„+N  aj  -  (N  ayi(}i^  +  N  aj), 

}j.  +  N 

4^  =  NaJ(p„+  Na_). 


in.  B.  3.  Strong-Property-Fluctuation  Theory  (SPFT) 

An  alternative  to  the  various  forms  of  the  pronounced  Maxwell  Garnett  model  is 
Strong-Property -Fluctuation  Theory  (SPFT).  This  theory  is  based  on  variations  of  the 
Maxwell  equations  and  position  averaged  values  of  individual  particle  polarizabihties. 
(SPFT)  has  been  formulated  for  homogenizing  chiral  composites  [3].  The  argument  begins 
with  the  Drude-Bom-Fedov  constitutive  relations  of  a  nonhomogeneous  chiral  medium 

D(r)  =  e  (r)[E(r)  +  fi  (r)V  x  E(r) , 

B(r)  =/z  ir){Hir)+p  (r)V  xH(r), 


where 

e  (r)  =  permittivity  scalar  ^ 

(r)  =  permeability  scalar  I  All  three  being  imphcit  functions  of 
P  (r)  =  chirality  pseudoscalar  )  (O  (circular  frequency) 

Rewritten,  the  constitutive  relations  in  compact  matrix  notation  are  given  by: 

C(r)  =  (r)  F(r) 


with 

C(r)=  fD(r)V 
l,B(r)J 


fO  i^ 
0)  l^-i  Oj 


j/(r)=  7+(r)  7-(r) 


(P  (r)  iCQUij)  ^ 

I  k(rf  I 

\-i  Q)e(r)  P{r)  \  ,  F(r)  =  fE(r)^ 
I  k{rf  )  lH(r)j 
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where 

k  (r)^QrJe  (r)  [i  (r)  and  7  +(r)  =  k  (r)/[  1-  yt  (r)  ji  (r)] 

7-(r)  =/:  (r)/[l+^  (r)  p{r)\ 

I^e  goal  of  this  theory  was  to  generalize  the  strong-permittivity-fluctuation  theory  of 
^electoc  mixtures  to  chiral-chiral  composites,  a  composite  material  in  which  the  particulate, 
in  addition  to  the  matrix  media,  are  chiral.  A  detailed  explanation  of  this  theory  can  be 
found  in  reference  [3]. 

Up  to  this  point,  the  constimtive  relations  presented  have  only  applied  to 
variations  homogeneous  or  inhomogeneous  forms  of  isotropic  chiral  media.  The  following 
constitutive  relations  will  form  the  basis  for  the  behavior  of  linear  polarized  waves  in 
anisotropic  chiral  composites. 

in.  C.  Electromagnetic  Fields  in  a  Composite  Anisotropic  Chiral  Medium 

A  homogeneous  composite  anisotropic  chiral  material  is  described  by  the  following  tensor 
constitutive  relations 


where 

£ 


D  =£  •E-hi^B 

H  =  E  +  •  B 

-ig 

o1 

1  ig 

e 

0 1 ,  is  the  permittivity  tensor 

Lo 

0 

-ik 

o1 

1  ik 

0  I ,  is  the  permeability  tensor 

L  0 

0 

and  is  the  chirality  admittance  [20]. 


Fenite  is  any  of  several  magnetic  substances  that  consist  essentially  of  an  iron 
oxide  combined  with  one  of  more  metals  (as  manganese,  nickel,  or  zinc),  have  high 
magnetic  permeability  and  high  electrical  resistivity.  Another  proposed  form  of  composite 
is  a  combination  of  chiral  and  ferrite  media. 

Constitutive  Relations  for  Nonreciprocal  Composite  Chiral-Ferrite  Media 

The  constitutive  relations  for  chiral  media  are  [19] 


D  =€E  +^H, 
B  =//H  +^E, 
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+  ia^.j 
C  —  ^cr2  ^^ci2 

where  for  lossless  chiral  media 

^crj  =  ^cr2  =  ^  rcciprocal  case), 

^  cil~  ~^ci2 

^  0,  ^  0  (for  the  nonreciprocal  case), 


The  constitutive  relations  for  ferrite  media  are  [19] 

D  =  £E, 

B=/lH, 

where 

r  ^  iK  oi 

ti  =  \-ik^  0  I , 

L  0  0  /to  J 

/I  =  /j.‘  - /i" ,  k„=  k„‘  -ikj’,  and  the  dc  magnetic  field  B„  is  in  the  z  direction. 

Thus,  the  composite  medium  is  hypothesized  to  have  the  following  constitutive  relations: 

D  =  eE  + 

B=fiH  +  CE. 

in.  D.  Wave  Propagation  and  Scattering  Characteristics  of  Chiral 
Materials 

The  wave  propagation  and  scattering  characteristics  of  chiral  materials  differ  from 
those  of  G/E  composites  (nonchiral  composites)  because  the  conversion  of  linear  polarized 
waves  to  alternating  eUiptical  RCP  and  LCP  ,  right  circularly  polarized  and  left  circularly 
polarized  fields,  respectively.  The  difference  depends  on  the  surface  induced  current  and 
charge  densities  at  the  boundary  of  chiral  materials.  The  RCP  and  LCP  fields  produced  by 
the  incident  linearly  polarized  waves,  within  the  spherically  shaped  chiral  inclusion,  cause 
differences  in  velocity.  If  multiple  interfaces  between  the  chiral  inclusion  and  the  achiral 
matrix  are  incorporated  into  the  design  of  the  chiral  material,  multiple  scattering  occurs  due 
to  alternating  conversion  patterns,  i.e.  RCP  to  LCP  fields  and  LCP  to  RCP  fields.  Then, 
as  the  system  experiences  multiple  scattering,  the  decrease  in  the  effective  transmission  of 
em  waves  will  increase  absorption  [16]. 
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in.  D.  1.  The  Small  Spheres  Assumption 

The  use  of  spherical  geometry  when  referring  to  the  individual  inclusion  particles, 
whether  hehcal  or  omega  (Q)  shaped,  is  assumed  for  uniform  size  distribution.  This 
assumption  reduces  comphcations  due  to  inclusion  size  and  shape.  Further  studies  have 
shown  that  the  absorption  patterns  due  to  cross  section  are  larger  for  spheres  than  those  for 
spheroids.  For  this  reason,  the  term  chiral  sphere  is  used  when  addressing  the  chiral 
inclusions  embedded  in  an  achiral  matrix  [16]. 

III.  D.  2.  Transmission  Line  Theory  For  Chiral  Composites 

Initially,  the  transmission  line  theory  was  used  to  develop  an  understanding  of  how 
an  electromagnetic  wave,  in  the  form  of  a  plane  wave,  propagates  through  a  composite 
plate  composed  of  a  series  of  uniformly  constructed  lamina.  In  this  paper,  the  developed 
theory  was  termed,  “The  Foundational  Theory”.  Assuming  the  same  criteria  on  a 
macroscopic  basis,  a  T  matrix  can  be  formulated  for  spherical  chiral  inclusions  embedded 
within  an  achiral  matrix.  The  final  result  of  a  T  matrix,  presented  below,  for  a  single  chiral 
sphere  in  nonchiral  (or  achiral  )  dielectric  medium.  The  boundary  conditions  on  the 
electromagnetic  fields  ( E-  and  H-fields)  inside  and  outside  the  sphere  (assuming  continuity 
at  the  interfaces)  which  are  expanded  through  vector  spherical  functions  such  as  M(3)  and 
N(3),  in  the  T  matrix  and  LCP  and  RCP  fields  in  a  chiral  composite  (which  is  analogous  to 
the  A;  and  B-  coefficients,  where  i=l,  2,  3,.  .  .  ,  n  number  of  lamina  or  vectors  for  the 
filament-current  phase  correction  model )  [14], [16]. 

Thus,  the  analytical  expressions  for  the  effective  dielectric  constant  can  be  obtained  from  a 
software  package  called  MACSYMA  .  Once  the  analytic  expressions  for  the  lowest  order 
T-matrix  elements  are  known  in  long-wavelength,  results  can  be  used  in  the  Maxwell 
Garnett  and  Bruggeman  models  to  calculate  initial  estimates  to  the  roots  of  dispersion 
equations  at  higher  frequencies  and  concentrations  of  spherical  particles  [16]. 
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where  the  elements  of  the  T  matrix  for  a  chiral  sphere  in  a  dielectric  host  medium  are  given 
by 


CT'  \on=  -  [W„(L)A„  (R)  W„  (R)  A„  (L)]/[  W„  (L)V„  (R)  W„  (R)  V„  (L)] 
cr'")een  =  "  [B„  (L)V„(R)  -h  B„  (R)  V„  (L)]/[  W„  (L)V„  (R)  +  W„  (R)  V„  (L)] 
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( T'  ■)„.  =  -  [W„(L)B.  (R)  -  W.  (R)  B,  (L)]/[  W.  (L)V.  (R)  +  W.  (R)  V.  (L)] 

where 

W„  (J)  =  m(p„(mJx)(C  (x)  -  (p„’(mJx)  ^  „(x) 

V„  (J)  =  9n(mJx)(C  n’  (x)  -  m(p„’(mJx)  ^ „(x) 

A„  (J)  =  (p„(mJx)  (p/(x)  -  m(p„’(mJx)  %  (x) 

B„  (J)  =  m(p„’(mJx)  (p„’(x)  -  (p„’(mJx)  (p„  (x) 

J  is  either  L  or  R,  and  x  is  k’a  (where  the  prime  denotes  differentiation  with  respect  to  the 
argument).  (p„  (z)  =  z  j  „  (z),  ^  „(z)  =  z  h  „  (z),  and  j  „  and  h  „  are  the  spherical  Bessel  and 
Hankel  functions,  respectively.  The  parameters  m^  =  k^/k’,  mj^  =  1%/k’,  and  m  = 

2mLmjj/(mL  +  mjj)  [16]. 

in.  D.  4.  Propagation  in  Anisotropic  Chiral  Medium 

The  general  bianisotropic  medium  is  characterized  by  the  constitutive  equations  [4] 

D  =  e*E  +  ^*H, 

B  =  [i*H  +  C  -E 


Assuming  lossless  reciprocal  medium. 

Now,  assuming  isotropic  permittivity  and  permeability,  the  constitutive  relations  for  the 
general  anisotropic  chiral  medium  follows 

D  =  8E  -jV(ep)K^*  H, 

B  =  pH  +  jA/(Ep)K^«  E 

the  relative  chirality  dyadic  being  real.  Consider  a  special  case  where  a  chiral  medium  is 
made  of  metal  helices  with  different  handedness  aligned  with  x  and  y  directions.  The 
chirality  dyadic  can  be  taken  as 


\  =  S  (‘kUx  -  UyUy). 

The  wave  numbers  are 


K  =  k _ Vi-k!^ 

Vi  ±  sin^G  cos2(p 
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In  the  present  case,  there  are  two  planes  defined  by  cp  =  ±  7t/4.  The  wave  numbers  for  the 
right-handed  and  left-handed  ellipsoid  surfaces  illustrated  in  Figure  1  are  k-  (0,  tp)  and  k-i- 
9)’  respectively.  The  values  of  the  main  axis  of  the  ellipsoids  are 

K,,=kVriH±)i^i,  K^,=kVri±Kj,  K^=kV[i-Kj. 

Li±s-I  Li(-/+)k^J 


Figure  1-3.  The  Elliptic  Wave  Number  Surfaces  of 
k-(0,(p)  and  k+(0,(p) 
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in.  E.  Constitutive  Relations  for  Omega  Shaped  Inclusions 

While  advancing  the  concept  of  helical  shaped  inclusions  in  chiral  composites, 
omega  shaped  inclusions  demonstrated  properties  which  enhanced  absorption  and/  or 

reduced  reflection  from  planar  lossy  slabs.  Metal  particles  shaped  like  an  Q  could  reduce 
the  reflection  of  surfaces  when  illuminated  by  linearly  polarized  and  mipolarized 
electromagnetic  waves.  The  stems  of  the  particles  would  align  with  the  x-axis  with  the  loop 
in  the  (x-y)  plane.  The  medium  which  consists  of  these  omega  particles  are  termed  uniaxial 
bi-anisotropic,  with  constitutive  relations  given  by  [17] 


D  =£-E  +  iK^J(  •  H 

B  .  H  +  iK^/(  ‘  E 


Here 

The  uniaxial  dyadics  with  transverse  (t)  and  normal  components  (n)  are: 


where  z„  stands  for  the  unit  vector  along  the  z-axis  (normal  to  the  interfaces)  and  the 
transverse  unit  dyadic  and  the  90°  rotator  in  the  (x-y)  plane  are: 


I,  =  +  y„y„  and  J  =z„  xl,  =  y^„  - 

Additional  coupling  provided  by  the  omega  particles  is  measured  by  the  coupling  parameter 
K,  given  by  the  following  equation: 


K 

2 

when  the  reflection  coefficient  equals  zero.  The  respective  reflection  and  transmission 
coefficients  for  a  nonsymmetric  slab  in  air  can  be  found  using  transmission-line  theory  and 
are  expressed  as  the  following  [17]: 


R  =  _ (Y^  -i-y-irT,  -Y+)fl-exv(-2jpdJl _ 

(Y^  +  Y+)(Y„  +  Y-)  -  (Y„  -  Y-)(Y„  -  7+)  exp  {-2m 


T  =  _ 2^Y+  +Y.)Y  exx^(-iBd) 

(Y„  +  Y+)(Y^  -t-  Y-)  -  (Y„  -  Y-)(Y^  -  T+)  exp  (-2jpd) 

where  d  is  the  slab  thickness  and  Y^  symbolizes  the  free  space  wave  admittance. 
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The  polarization  of  transverse  electric  and  magnetic  fields  are  best  characterized  by  wave 
impedances  or  admittances  for  plane  linearly  polarized  electromagnetic  waves  propagating 

in  uniaxial  Q  media.  The  admittances  for  TM-  and  TE-polarized  fields  are: 

Ti:=V  /£■£■  )  1  xpsJ(l-_Jcf_  -^)(-l+)j  K__l 

L  e//,  Vfe/t, )  J 

ko'ej^. 

Y±  =  V/:44  _Af _  -K^  )  (-!+}]  K__l 


respectively,  where  the  upper  sign  corresponds  to  the  waves  in  the  positive  z-direction  and 
the  lower  sign  to  the  waves  traveling  in  the  negative  z-direction.  Here  k^  is  the  transverse 

part  of  the  propagation  factor  and  -  0>J(e^iiJ. 

The  propagation  factors  of  the  TM-  and  TE-polarization  waves,  /?  ,  included  in  the 
reflection  and  transmission  coefficients  are: 

P  =  ±Vf£,_(k/ej2,  -kf  )-  kfK^J  and  P=  -k^  )-  kfK^  ] 

respectively  [17]. 

Therefore,  there  are  three  material  parameters  that  control  the  reflection  and  transmissions 
properties  of  uniaxial  omega  bianisotropic  composite  slabs:  the  dielectric  permittivity, 
magnetic  permeability  and  the  coupling  parameter.  Thus,  the  adjustment  of  the  magneto¬ 
electric  interaction,  K  ,  the  input  impedance  and  the  free-space  impedance  can  be  matched 
to  achieve  low  levels  of  reflection  over  wide  frequency  bands  [17]. 
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IV.  other  Composite  Types 


IV.  A.  Quasioptic  Composites 

The  Faraday-Kerr  effect  is  the  scientific  explanation  of  how  certain  material 
properties  change  plane  polarized  light  at  optical  frequencies  (Figure  1-4).  By  using 
composites  made  of  rare-earth  materials  and  ferrites,  similar  effects  can  be  simulated  at 
quasioptical  frequencies,  otherwise  known  as  microwave  and  rnillimetCT  wave  frequencies. 
Such  composite  types  are  useful  for  changing  the  polarization  of  incident  microwave 


Figure  1-4.  The  Faraday  Effect 
describes  the  rotation  and  the  Kerr  Effect 
explains  the  directionality.  However  the 
Kerr  Effect  can  be  described  by  three 
different  subeffects: 


a)  Polar  effect  -  Magnetization  occurs 
normal  to  the  surface  of  reflection. 

b)  Longitudinal  effect  -  Magnetization 
occurs  parallel  to  the  surface  of  reflection. 

c)  Transverse  effect  -  Magnetization 
occurs  perpendicular  to  the  plane  of 
incidence. 


One  of  the  major  concerns  in  the  manufacture  of  composites  is  low  cost,  high 
quality.  The  design  of  modem  electronic/communication  components  requires  low-loss, 
high-permittivity  dielectrics.  For  example,  consider  a  two-phase  mixture  of  polystyrene,  a 
high-quality  dielectric  with  low  permittivity  and  a  dielectric  constant  of  2.1;  and  barium 
titanate,  an  expensive  high-permittivity  material,  with  a  dielectric  constant  greater  than 
1,000.  With  the  appropriate  volume  fractions,  the  composite  formed  would  possess  the 
low-cost,  low-permittivity  properties  of  polystyrene  and  the  high-permittivity  properties  of 
barium  titanate.  Thus,  the  resulting  composite  would  have  high  permittivity  with  respect  to 
polystyrene  and  low  cost  with  respect  to  barium  titanate  [1]. 


energy  [1]. 


References: 
http:Hphsg.uni- 
bielefeld.del  helium/ carsil.htm 

http:llbessrc.msd.anl.govlllid.htm 


45 


IV.  B.  Dielectric  Composites 
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CONCLUSION 


In  conclusion,  graphite/epoxy  composites  currently  make  the  greatest  contribution 
to  industrial  and  military  systems.  Chiral  composites  have  the  potential  to  make  substantial 
contributions  to  future  applications  due  to  the  additional  degree  of  freedom  supplied  by  the 
chirality  parameter.  Finally,  quasioptic  and  dielectric  composites  have  electromagnetic 
properties  that  produce  the  predictable  behavior  needed  in  electromagnetic  radio  frequency 
responses.  A  separate  research  study  is  suggested  for  the  latter  composites  due  to  their 
potential  economic  and  structural  improvements.  Additional  research  topics  include:  1) 
Techniques  used  to  Measure  Parameters  of  Advanced  Composite  Materials;  and  2) 
Evaluation  of  Advanced  Composite  Materials  as  a  Function  of  Frequency  and  Temperature. 
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MISSION 

OF 

ROMELABORA  TORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


